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Abstract. We present the first maps of liigh-excitation CO J— 6-5 and 7-6 
and isotopologue lines over ~ 2' — 5' regions at 10" resolution toward low-mass 
protostars to probe the origin of the warm gas in their surroundings. The data 
were obtained using the CHAMP+ 650/850 GHz heterodyne array receiver on 
APEX0 Surprisingly strong quiescent extended narrow-line high-J ^^CO 6-5 
and 7-6 emission is seen toward all protostars, suggesting that heating by UV 
photons along the outflow cavity dominates the emission. At the source position 
itself, passive heating of the collapsing inner envelope by the luminosity of the 
source also contributes. The UV photons are generally not energetic enough to 
dissociate CO since the [C I] 2-1 emission, also probed by our data, is weak 
except at the bow shock at the tip of the outflow. The extended UV radiation is 
produced by the star-disk boundary layer as well as the jet- and bow-shocks, and 
will also affect the chemistry of species such as H2O and HCN around the outflow 
axis. Shock-heated warm gas characterized by broad CO line profiles is seen only 
toward the more massive Class outflows. Outflow temperatures, estimated 
from the CO 6-5/3-2 line wing ratios, are ~100 K. These data illustrate the 
importance of getting spatial information to characterize the physical processes 
in YSO surroundings. Such information will be important to interpret future 
Herschel and ALMA data. 



1. Introduction 

During the formation of low-mass stars, material moves from the cloud core to 
the collapsing envelope and disk onto the growing star. At the same time, part 
of the envelope is dispersed by the jets and winds from the protostar, limiting its 
growth. Characterizing these different physical components of low-mass young 



^This publication is based on data acquired with the Atacama Pathfinder Experiment (APEX). 
APEX is a collaboration between the Max-Planck-Institut fiir Radioastronomie, the European 
Southern Observatory, and the Onsala Space Observatory. 
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stellar objects (YSOs) (envelope, disk, outflow, cloud) is thus crucial for under- 
standing the physical evolution of the earliest stages of star formation: this is 
the phase in which the final mass of the star and disk are determined and the 
conditions for planet formation are set. Because of the high extinction (Ay > 50 
mag), these sources are best studied at (sub)millimeter and far-infrared wave- 
lengths. Most studies to date have focussed on the cold gas and dust around 
protostars, in particular using the dust co ntinuum and low- J l i nes of CO and its 
isotopologues ^^CO, C^^O and/or C^' ^O (jShirlev et ani2000l : Uo^sensen et al. ' 



[2002 . .20051 : iHatchell k Dunham \\200^ ). In contrast, the warm gas is much more 
diagnostic of the energetic processes that shape these deeply embedded sources. 

The most direct probes of warm (50-200 K) gas are the high-J (J > 4) 
lines of CO (for example, Eu = 115 K for J=6). Complex molecules such 
as H2CO and CH3OH also have high excitation lines in the 230 and 345 GHz 
atmospheric windows which have been detected toward low-mass protostars (e.g. 
van Dish oeck et al. I|l995| : iBlake et ain il995l: ICeccare Ui et al. f2000l:|Schoier et al. 



2002; Ma ret et al. I |200j; iJgrgensen I 2004; Bottinelli_etalJ l2004i . l2007l ). but 
their strong abundance variations complicate their use as tracers of the phys- 
ical structure. Observations of high-excitation lines of CO can be carried out 
from the ground in the atmospheric windows at 650 and 850 GHz, but they 
require excellent weather conditions. Such pioneering data were obtained more 
than a decade ago at single positions with the Caltech Submillimeter Obser- 
vatory (CSO) using its excellent receiver suite ([HoEcrliciidc et al. , 19 98) and 
with the James Clerk Maxwell Telescope (JCMT) (|Schuster etlDfiggsl . Il995l l. 
Around the same time, the Infrared Space Observatory (ISO) provided spec- 
trally unresolved data on even higher- J l i nes in a large ^ 80^^ beam (e.g. 



Ceccarelli et"anil998l: IG lannini et al. II1999I : iNisini et aljl999l : ICiannini etaL 



2001; Nisini et al. 1 20021 ) . The lack of spatial and/or spectral information on the 



warm gas prevented an in-depth analysis of these data, however. 

As a result, the location and heating mecha nisms of warrn dense gas near 

low-mass protostars are s till strongly debated (e.g.. lNisini et al.l2000l : [Ceccarelli et al. 



I999I : iMaret et al. Il200"3 ). Proposed options include (i) passive heating of the 
inner ~100 AU region of the collapsing envelope by the protostellar luminosity; 
(ii) active heating in shocks created by the interaction of jets and winds from 
the protostar with the envelope and cloud out to large distances; (iii) heating 
by UV photons escaping through outflo w cavities and scatt ered back into the 
envelope on a few thousand AU scales (jSpaans et al. Ill995l): and (iv ) a form- 
ing protoplanetary disk heated by accretion shocks ( Ceccarelli et al."l[2002. ). In 
scenarios (i) and (iv), the high-J CO emission should be spatially unresolved 
and centered on the source. In scenario (ii), extended broad emission lines are 
expected to be seen out to the tip of the outflow. In scenario (iii), spatially ex- 
tended narrow high-J lines are expected. APEX-CHAMP"^, with its heterodyne 
spectral resolution, 7-10" beam (~ 1000 AU at 120 pc) and mapping capabili- 
ties, can directly test the various scenarios. Located at the Chajnantor plateau 
in Northern Chile with excellent atmospheric transparency, it allows much more 
routine observations of these high-J lines than possible from Mauna Kea. We 
report here the initial results of a survey of high— J CO lines toward a set of 
low-mass protostars (van Kempen et al. 2009a, b). 
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2. Observations 



The CHAMP+ instrument, developed jointly by the MPIfR and SRON Gronin- 
gen, is the only submillimeter array receiver in the world able to simultaneously 
observe molecular line emission in the 650 and 850 GHz atmospheric windows on 
arcminute spatial scales with 7-9" pixels ( Kasemann et al.1l2006l : iGiisten et al. I 
20081 ). It is the successor of the CHAMP array at 460 GHz deployed on the CSO 



until 2003. CHAMP"*" has 14 pixels (7 in each frequency window) arranged in a 
hexagon of 6 pixels around 1 central pixel. To obtain fully sampled maps, the 
array was either moved in a small hexagonal pattern, or, more efficiently, in an 
on-the-fly (OTF) mode. The hexa-pattern covers a region of about 30" x 30", 
whereas small OTF maps are 40" x 40", or occasionally larger. 

The observations were obtained during several observing runs between June 
2007 and November 2008 in three different line settings, with the combination of 
lines chosen to match required integration times: (1) ^^CO J=6-5 and 7-6; (2) 
^^CO 6-5 and [C I] 2-1; and (3) C^^O 6-5 and ^^CO 8-7. For the latter setting 
(so far done only for HH 46), a stare mode was used to increase the S/N on the 
central pixel. A position switch of 900" or larger was used for all settings, except 
for the stare setting (3), which used a beam-switching of 90". Complementary 
lower- J lines of CO, HCO"*" and their isotopologues have been obtained with the 
facility APEX-1 and APEX-2a receivers. 

During the observations in 2007, the backend consisted of two Fast Fourier 
Transform Spectrometer (FFTS) units serving the central pixel (resolution up 
to 0.12 MHz), and 12 MPI- Auto-Correlator Spectrometer (MACS) units (res- 
olution 1 MHz) connected to the other pixels. Since July 2008, all pixels are 
connected to FFTS units. Main beam efficiencies, derived using observations 
on planets, are 0.56 at 650 GHz and 0.43 for 850 GHz. Typical single sideband 
system temperatures are 700 K and 2100 K, respectively. Pointing was checked 
on various planets and sources and was found to be within 3". Calibration is 
estimated to have an uncertainty of 30% for both frequencies. Integration times 
were such that the typical rms in a 0.7 km s^^ velocity bin at 650 GHz in setting 
(1) is 0.2-0.4 K, in setting (2) 0.2 K, and setting (3) 0.1 K, respectively. At the 
map edges, noise levels are often higher due to the shape of the CHAMP"*" array. 

For comparison, the Herschel Space Observatory will allow observations 
of far-IR CO lines with the PACS instrument at spatial resolutions similar to 
these APEX data (~10"). At its longer wavelengths (~500 /um), the beam of 
Herschel is comparable to, or smaller than, the field of view of CHAMP"*" (~ 40"). 
Thus the CHAMP+ data provide important complementary information on the 
distribution of warm gas within the Herschel-HIFI beams. 

The sample studied so far consists of a dozen well-known and well-studied 
embedded protostars that can be observed from the southern sky. All sources 



have been studied in previous surveys of emb edded YSOs (e.g.. lJ0rgensen et al. 



20021 : iJoreensen I l2004l : iGroppi et al7ll2007l ^. Tabled] gives the parameters o:' 



ea ch source and its pr o perties . Initial CHAM P"*" results have been publi shed 
bv Ivan Kempen et al. I (|2009al ) for HH 46 and Ivan Kempen et al. I ( 2009bl ) for 
a larger sample. The data on NGC 1333 IRAS 4A/B (Yildiz et al., in prep.) 
and the Serpens core containing SMMl, SMM3 and SMM4 (Kristensen et al., in 
prep.) are being analyzed. Results using the single-pixel 460/850 GHz FLASH 



4 



van Dishoeck, van Kempen <fe Giisten 




10 10 

Velocity (km/s) 



CO 7-6 (x2) 



[C I] 2-1 (x5) 




-10 10 20 

Velocity (km/s) 



Figure 1. Left: Spectra of various CO and isotopologue lines centered at 
the source position. Right: Spectra at an off position in the red outflow lobe 
at (-20", -20"). The quiescent gas has a FWHM of 1.5 km s'^ centered 
at Vlsr=5.3 km s~^ (based on van Kempen et al. 2009a, published with 
permission) . 



receiver on APEX can be found in van Kempen et ITI tOOd . \2Q09(\ ) (see also 
van Kempen 2008). 



3. HH 46 as an example 



3.1. The source 

The first source targetted with APEX-C HAMP+ was HH 46, located at the edge 
of an isolated Bok globule {D = 450 pc) ([Schwartz 19771 ). It is well-known for its 
spectacular outflow, observed at both visi ble and infrared wavelengths with the 



Hubble and Spitzer Space Telescope (e.g., Heathcqte et al. I Il996l : IStanke et al. 



l eLesco vt . , , . . 

119991 : iNoriega-Crespo et al. Il2003 : Ivelusamv et al. ^00^ Deep Ha observations 
have revealed bow shock s associated with th e HH 46 outflow up to a parsec away 
from the central source ( Stanke et al. 19991 ). Its blue-shifted lobe expands into 



a low density region outside the cloud, whereas the red-shifted lobe plows into 
the den se core. The internal driving source for the flow is HH 46 IRSl (L = 
16Lm) (iRavmon d et al. I|l994l: ISchwartz fc Greene l[200^ ). IChernin &: Masson 



(119911 ) and Olbcr get al. I (|l992l ~ mapped this region using low excitation CO 



lines showing that the molecular emission in the red-shifted outflow lobe is indeed 
much stronger than that in the blue-shifted one. 
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Table 1. Sample of sources observed with CHAMP"*" as of summer 2009 



Source 


RA 


Dec 


D 


Lhoi 


7b ol 


Class 


Ref. 




(J2000) 


(J2000) 


(pc) 


(Le) 


(K) 






NGC 1333 IRAS2 


03:28:55.2 


+31:14:35 


250 


12.7 


62 





2 


NGC 1333 IRAS4A/B 


03:29:11 


+31:13 


250 


8 


35 





3 


L1551 IRS5 


04:31:34.1 


+ 18:08:05.0 


160 


20 


75 


1 


2 


TMR 1 


04:39:13.7 


+25:53:21 


140 


3.1 


133 


1 


2 


HH 46 


08:25:43.8 


-51:00:35.6 


450 


16 


102 


1 


1 


Ced 110 IRS4 


11:06:47.0 


-77:22:32.4 


130 


0.8 


55 


1 


2 


BHR 71 


12:01:36.3 


-65:08:44 


200 


11 


60 





2,5 


IRAS 12496-7650 


12:53:17.2 


-77:07:10.6 


250 


24 


325 


1 


2 


Serpens core 


18:29:55 


+01:14 


250 


30/5 







4 


RCrA IRS7 


19:01:55 


-36:57:21 


170 
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References: (1) van Kempen et al. (2009a); (2) van Kempen et al. (2009b); (3) Yildiz 
et al. in prep.; (4) and (5) Kristensen et al. in prep. 



3.2. Results 

Figure [1] (left) shows the high-quality spectra taken at the position of HH 46 
IRSl. Emission is clearly detected for all lines, except for C^^O J=6-5 and 
^^CO 8-7. [C I] 2-1 emission is surprisingly weak at the central source, only 
^MB = 2.7 K. Line wings are clearly seen for the ^^CO lines but become less 
prominent relative to the narrow line component for the higher excitation lines. 
The same is seen at outflow positions away from the source (Fig. [H right) 

Figure [2] presents the ^^CO 6-5 spectra binned to square 10" x 10" pixels. In 
addition to the line wings characterizing the swept-up outflow gas, surprisingly 
strong extended narrow CO lines are seen along the direction of the outflow 
axis (see for example spectra at —30", 10" and —30", —30"). These data directly 
support scenario (iii), in which UV radiation heats the surrounding outflow walls. 

To test this scenario further, the outflow region was mapped more exten- 
sively in ^^CO 6-5 and [C I] 2-1 lines out to the tip of the outflow characterized 
by the bow shock around (—100", —60"). Figure [3] shows an overlay of the CO 
and [C I] integrated intensity maps on the Spitzer image. The CO emission (both 
narrow and broad) clearly disappears toward the location of the bow shock. Al- 
though the dynamic range of the [C I] map is not large, narrow [C I] 2-1 line 
emission is somewhat enhanced around the bow shock position. This suggests 
that along most of the outflow axis the UV photons are energetic enough (>few 
eV) to heat the gas, but not to photodissociate CO. In contrast, at the bow- 
shock position, photons must be present with high enough energies to dissociate 
CO (>11 eV, <1100 A) and produce [C 11 emission both ahead and just behind 
the shock. According to the models bv iNeufeld Sz Dalgarno I ( 19891 ). J— shocks 



with velocities larger than 90 km s ^ produce CO dissociating photons, mostly 
through the two-photon decay of metastable He(2 ^S) atoms. At the bow shock. 



the s hock velocity is measured to be 220 km s ^ from atomic lines ( Fernandes 



200d ). more than sufficient to produce these high energy photons. Because of 



the high densities in the core, they can travel only a short path before being 
absorbed by dust, H2 or CO. This situation is reminescent of the IC 443 su- 
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Figure 2. Spectral map of the ^^CO J=6-5 lines obtained with APEX- 
CHAMP+. The thick line indicates the outflow direction. Note the narrow 
CO line profiles offset from the outflow walls, in addition to the broad out- 
flowing gas (see blow-up of spectrum at —50", —60") (based on van Kempen 
et al. 2009a). 



pernova remnant, where quiescent [C I] emission was also detected ahead of the 
shoc k and interpreted t o arise from photodissociation of CO in the pre-shocked 



gas (|Keene et al. Ill996l ). 



The CHAMP+ data thus prov ide support fo r scen ario (iii), but with an 
extension of the original model by ISpaans et al. I ( 19951 ). First, the UV pho- 
tons originate not only from the sourc e itself but are also p roduced further out 
by the bow- and jet-shocks. Second, ISpaans et al. I ( 19951 ) suggested that the 
UV photons from the star-disk accretion boundar y layer reach large dist ances 
by scattering on dust grain s in the outflow cone. iBruderer et al. I (|2009l ) and 



van Kempen et al. I ( 2009al ) note that direct illumination of the outflow walls 
may be more effective if the outflow walls have a parabolic rather than a linear 
straight shape (see Figured] for illustration). 

To quantify the contribution to the emission at the source position from 
scenario (i), a pas sively heated env e lope model has been constructed follow- 
ing the method of |j0rgensen et al~ ( 20021 ). The envelope is assumed to have 
a power-law density structure, with the temperature determined by a dust ra- 
diative transfer calculation with the luminosity of the source as input. The 
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Figur e 3. Spitzer three co l or (3. 6 (blue), 4.5 (green) and 8 /im (red)) image 
froni iNoriega-Crespo et aTl ()2004f ) with the contours of integrated CO J=6- 
5 emission (blue/dark) observed with APEX-CHAMP+ overlaid. The CO 
contours are in steps of 5 K km s~^. The [C I] 2-1 emission (green/light) is 
detected weakly on source but peaks further down at the bow shock at the 
tip of the outflow where the UV photons produced in the fast bow shock are 
hard enough to dissociate CO. Analysis of the line profiles shows that the 
emission consists both of accelerated swept-up gas along the outflow as well 
as quiescent, photon-heated gas surrounding the outflow cavity walls (see also 
Fig. 2) (reproduced with permission from van Kempen et al. 2009a). 



power-law exponent, mass and extent of the envelope are determined by a best 
fit to the spectral energy distribution and dust continuum maps. For HH 46, 
a LABOCA 850 /xm map over the central 5' x 5' map has been obtained for 
this purpose. Subsequently, the CO excitation and line emission are calculated 
for the resulting best-fit temperature and density envelope structure us i ng th e 
Monte Carlo radiative transfer package of iHogerheijde Sz van der Tak I (|2000l '). 
assuming that the gas temperature equals the dust temperature. The latter 
assumption should hold for the inner envelope where densities are high enough 
(> 10^ cm~'^) for gas and dust to be thermally coupled. 

This quantitative analysis shows that the envelope emission can reproduce 
some of the observed high- J CO emission shown Figure [1] (left), but generally 
underproduces the ^^CO and ^^CO 6-5 lines by nearly a factor of 3. Increasing 
the CO abundance is not possible since then the C^^O 3-2 and 6-5 lines would be 
overproduced. Thus, an additional optically thin but warm component needs to 
be added to the model to explain the data. It needs to be relatively unobscured 
since it adds emission to that of the optically thick ^^CO lines from the inner 
envelope. Photon-heated gas in the (outer) envelope fulfills these criteria. The 
weak observed [C I] emission at the source position can be accounted for by a low 
atomic C abundance (-^few xlO"''' with respect to hydrogen) that is maintained 
by cosmic ray induced UV photodissocation of CO inside the dense envelope. 
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Figure 4. Cartoon of the HH 46 outflow on scales up to the bow shock at 
^60,000 AU from the source, illustrating the photon heating of the cavity 
walls (either directly or scattered by dust) responsible for the bulk of the qui- 
escent high- J CO emission. UV photons are created at the star-disk accretion 
boundary layer as well as in the bow- and jet shocks. The inserts illustrate 
possible geometries at the base of the outflow, with the parabolic shape allow- 
ing more direct illumination of the outflow walls (reproduced with permission 
from van Kempen et al. 2009a). 



The current data do not provide constraints on scenario (iv), i.e., whether 
some of the hot gas could be due to the accretion shock onto the disk in the 
embedded phase. To answer this question, submillimeter interferometry in the 
CO J=6-5 and 7-6 lines with ALMA is needed. High spatial and spectral 
resolution mid-infrared CO v=l-0 data at 4.6 //m can also provide insight into 
this question. 

In summary, the quantitative analysis bv lvan Kempen et al. ( 2009al ) shows 
that the warm gas traced by the high- J CO arises from three different compo- 
nents: (i) the quiescent inner warm envelope, accounting for roughly 30% of 
the ^■^CO and ^^CO on-source intensities; (ii) the outflowing gas, accounting for 
the extended red and blue wings; and (iii) UV-heated quiescent gas (T > 100 
K), accounting (surprisingly!) for the hulk of the extended narrow 6-5 and 7-6 
emission. 
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Source: N1333IR4AB Line: C0(6-5) Freq: 691.473062 GHz Beam: 9.58 x 9.58 PA 0° 




Absolute coordinates Velocity (km.s ) 



Figure 5. Preliminary CO J=6-5 map of the NGC 1333 IRAS 4A/B region. 
The top map shows a channel image at the central velocity whereas the lower 
map shows an integrated intensity map. Two representative spectra are shown 
as well (Yildiz et al., in prep.). 



Other sources 



Data fo r several of the other sources in Table[T]are presented in lvan Kempen et al 



( 2009bl ) and confirm the results found for HH 46. Warm g SiS^ clS traced by ^^CO 



6-5 and 7-6, is present in all observed protostars at the central position. For 
Ced 110 IRS 4, a quantitative model similar to that for HH 46 shows that pas- 
sive heating of the envelope is again insufficient to explain the observed ^^CO 
6-5, 7-6 and ^^CO 6-5 lines, requiring heating of the envelope by UV photons 
even at the (0,0) position. 

Mapping shows that photon heating of the cavity walls takes place on arcmin 
scales for the outflows of several other Class and Class I protostars, as seen 
at positions off'-source of BHR 71, NGC 1333 IRS2, and L 1551 IRS5. Figure 
shows a preliminary map of the NGC 1333 IRAS4A/B region, with both outflow 
and narrow emission present along the outflow axis of the IRAS4A flow (Yildiz 
et al., in prep.). As for HH 46, the necessary UV photons are likely created by 
internal jet shocks and the bow shock where the jet interacts with the ambient 
medium, in addition to the disk-star accretion boundary layer. The distribution 
of the quiescent CO 6-5 and 7-6 emission seen toward BHR 71, where the narrow 
emission is stronger at larger distances from the source, confirms the hypothesis 
that UV photons necessary for heating can originate from both mechanisms. 
Although difficult to disentangle in cases where the outfiow is not in the plane 
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of the sky, the observations suggest that photon heating is a common phenomena 
along all outflows. The lack of [C I] 2-1 emission in the outflows indicates that 
there is limited production of energetic CO dissociating photons in the shocks. 

Shocked broad ^^CO 6-5 and 7-6 lines are only seen prominently in the 
flows of more massive sources (NGC 1333 IRAS2 and IRAS4A/B, BHR 71, 
HH 46, RCrA IRS7 and the Serpens sources), while lower mass flows do not 
show signiflcant^^CO 6 -5 and 7-6 line wings. iBachiller Tafalla I (Il999l ) and 
Arce Sz Sargent ( 20061 ) show that the outflows of Class I sources are more 



evolved and are driven with much l ess energy, producing significantly weaker 
shocks and less swept-up gas (see also iBontemps et al. Ill996l : iHogerheiide et alT 



19981 ). This is also reflected in the decreasing maximum velocity of ^^CO 6-5 



with lower outflow force. From the CO 6-5/3-2 line ratios, kinetic tempera- 
tures of ~100 K are found for the molecular gas in the flows studied here. Such 
te mperatures agree with expected conditions of bow-shock driven shell models 
bv lHatchell et al~ ( 19991 ). in which the momentum-conserving shells expand and 
loose kinetic energy to heat the swept-up molecular gas. Only the flows of L 1551 
IRS5 and IRAS 12496-7650 appear to be colder (<50 K), consistent with the 
'fossil', empty nature of these flows. 

The intensities at all positions toward RCrA IRS7 are an order of magnitude 
higher than can be produced by passive heating, and probably result from a 
significant PDR region near the source. Likely, RCrA itself heats the outside of 
the RCrA IRS7 envelope and cloud region. 



5. Concluding remarks 

Maps in high-J CO submillimeter lines clearly provide crucial information to 
determine the location and origin of warm dense gas near YSOs and the feedback 
that the young star has on its surroundings. Questions to be addressed with 
future larger samples combined with Herschel data include: is the UV-heated, 
quiescent component always dominant, or does this change from the deeply 
embedded Class phase to the less embedded Class 1 phase and eventually 
the disk-dominated phase? Does this heating limit the ability of the cloud to 
collapse further and thus help determine the final mass of the star? How hot is 
the outflowing molecular gas and how does this depend on evolutionary state? 

The processes identified here have strong implications not only for the phys- 
ical evolution, but also for the chemistry near protostars, especially the origin 
of the emission from water and complex organic molecules. For example, if UV 
radiation is indeed widespread around low-mass YSOs, it will readily dissociate 
H2O into OH and H, strongly affecting the predictions and interpretation of 
upcoming Herschel-HIFI data. Some indications of this process have already 
been found through enhanced narrow CN emissio n (a photoproduc t of HCN) 
along the outflow axis of the low- mass YSO L483 ( J0rgensen 20041 ) . In addi- 



tion, the UV photons can photoprocess and photodesorb ices such as CH3OH 
efficiently and thus enhance their gas phase abundances (Oberg et al. 2009a, b). 
Recent mapping of the CHAMP+ sources in HCN, CN and CH3OH with the 
JCMT-HARP and IRAM 30-m HERA arrays is in progress to test these chemical 
consequences (Kristensen et al., in prep). 
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Together, the quahty and richness of these data sets are a tribute to the 
instrument builders and illustrate how far and rapidly the field has developed 
since the first pioneering 230 GHz data bv lPhillips et al. I (|l973l ). 
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